Spectral flow in superconducting point contacts 
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' We find that multiple Andreev reflections mediating the transport in superconducting point con- 

tacts are strongly affected by a small amount of impurities in the area of the contact. We also argue 
that the model based on Zener transitions within independent conducting channels is not suitable 
C"j ' for kinetic processes in multichannel contacts. 
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INTRODUCTION 

Effect of impurities on transport in superconducting nanostructures is one of the key issues in the physics of 

| mesoscopic superconductors, among which point contacts are the simplest devices of interest. One can distinguish 
!h ' contacts with only a few transverse eigenmodes (or conducting channels) relevant for a given temperature (quantum 

contacts) and those where the number of channels N c h is relatively large (mesoscopic contacts). In quantum contacts, 
^ ■ the effect of impurities on both static and transport properties can be satisfactorily described by introducing certain 

1 t | transmission and reflection probability D + R = 1 for each conducting channel 0. Mesoscopic weak links with disorder 
3 ■ are sometimes also treated within this scheme. To this end, a quantum eigenstate problem for a charge carrier is to 

be solved for the particular realization of disorder, each state being considered as an individual channel with certain 
transmission probability. The problem of a multichannel weak link is then mapped into that for a set of tunnel 
contacts having an effective distribution of transmission probabilities. This scheme has been shown recently j2| to 
be equivalent to the semi-classical theory of superconductivity for calculating some static properties of multichannel 
O \ contacts. 

The individual channels found by solving the quantum-mechanical problem are independent in the sense that 
^vq [ transitions between them are absent under static conditions. One might now raise a question how useful is the model 
^ • of independent channels for transport properties of multichannel contacts when time-dependent perturbations induce 
transitions between the individual channels. This is of a crucial importance for further theoretical studies of transport 
in superconducting mesoscopic multichannel devices in presence of disorder. We address this problem using a specific 
| example of a ballistic point contact between two clean superconductors assuming that the thickness d of the insulating 
H ' layer and the size a of the orifice in it are much larger than atomic dimensions but shorter than the coherence length 
t-H \ £ and the mean free path £, i.e., a ~ d <C £ <C £. The contact thus has a very large number of conducting channels, 
^Voh ~ (pfo) 2 1. The spectrum of the midgap states localized within the coherence length £ near the point contact 
' was first found in 

S e± = T|A| cos(<£/2) (1) 

for a phase difference <p between the two superconductors. The spectrum is degenerate: all particles flying from one 
superconductor into another have the same energy independently of their momentum directions within the corre- 
] sponding hemispheres. Another example is a contact having a tunnel barrier with a transmission probability D. The 
• • ■ spectrum Q 
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e = ±|A|yi - Dsm 2 ((f>/2) (2) 

has a (mini)gap \AR|A| near the phase difference (f> = it, which separates the two branches of the spectrum. The 
mid-gap states have a profound effect on dynamics of the contacts. When the applied voltage is larger than the 
inelastic relaxation rate, eV 3> t~ x , a nonlinear regime sets in which originates from multiple Andreev reflections 
(MAR) [H IS IS IS E3 °f localized quasiparticles. For a contact without both impurities and a barrier, excitations 
from below the gap edge with e = — |A| are accelerated through the midgap states Eq. (JTJ and re- appear above the 
gap with e = |A| and vice versa, thus producing a highly noncquilibrium situation. On the contrary, if a contact has 
a barrier with a reflectivity R, the localized excitations cannot accelerate due to the minigap and retain their original 
states. It is only Zener tunnelling through the minigap that provides nonequilibrium current through the contact 

mm 



2 



The model of independent channels was applied to kinetic effects in multichannel point contacts in Ref. (£. The 
key ingredients of that model were (i) the MAR processes in separate channels with (ii) Zener transitions through a 
minigap between two spectral branches in Eq. J2J within each individual channel, and (iii) a standard distribution 
of reflection coefficients over the channels. Here we study the transport in mesoscopic multichannel contacts with 
impurities and show that the above model of individual channels consisting of processes (i)-(iii) is not suitable for 
dynamic phenomena when the number of channels is large, N c h 3> 1, because transitions between different channels 
are excited in time-dependent conditions. As a result, the system finds the optimal way to relax through channels 
with small reflection probabilities avoiding weak Zener tunneling. We further demonstrate that the concept of MAR 
is equivalent to the spectral flow appearing in a wide variety of quantum systems 0, EH, E3] • Using the quasiclassical 
theory of nonstationary superconductivity we show that the impurity scattering provides an effective momentum 
exchange between particles moving from the opposite sides of the contact and impedes the MAR acceleration. For 
low inelastic relaxation rate, the I-V curve displays a non-trivial behavior: At low voltages, impurity relaxation slows 
the acceleration down and results in a new linear current-voltage dependence with a large conductance. This differs 
from a square-root dependence obtained within the independent-channel model (i)— (iii). 



KINETIC EQUATION 

The density of states (DOS) for the midgap states in a classical point contact was calculated in Refs. 0, Q It is 
is(s) = i>(0)e~ k \ s \ with k — 2vp y|Ap — e 2 where s is the distance along the particle trajectory; is is normalized to 
the normal single-spin DOS at the Fermi level, Np. We use the standard kinetic equation approach [l5lH^| where the 
distribution function is split into two components /W and that are respectively odd and even in the variables 
e, p. They are coupled by the Boltzmann equation 

VVF 2Ts — + V — + VF dl K ) = J ■ (3) 

The second equation ensures that the distribution function is constant along the trajectory 0|. In Eq. J3J), the 
first term in the l.h.s. originates from the Doppler energy and x is the time derivative of the phase, and J is the 
impurity collision integral. A detailed discussion of the kinetic equations is given in Refs. 003. The standard 
technique of averaging over impurities applies since the number of impurities within the volume of the orifice is 
large. Indeed, the mean free time is r _1 ~ Npn- lmp \u\ 2 where \u\ ~ Up^ 3 is the Fourier transform of the impurity 
potential U. The number of impurities ni mp da 2 ~ (d/£)(Ep/U) 2 N c \ l can be very large even for U ~ Ep because of 
a macroscopic number of quantum channels in the orifice _/V c h 3> 1. The large number of impurities is the essential 
parameter in our theory. 

We integrate Eq. JSJ along the trajectory across the contact over distances longer than £. The first term in the 
l.h.s. of Eq. has a sharp maximum at s ~ d where is(s) is nearly constant. The third term in the l.h.s. disappears 
due to the decay of is. The collision integral contains |2( a combination of angular averages v \(is) — { v f^)\ • 
We write (is) = [(f) + + (f)_] /2 where (i/) ± = (27r) _1 J ±p >Q u(p, r) dfl p . The upper (lower) sign is for right-moving 
(left-moving) particles with p x > (p x < 0). We choose the direction of the x axis from the region with the phase 
Xi into that with X2- Since depends on the sign of p x only, being independent of the direction of p within each 
hemisphere p x > or p x < 0, 

J ± = -T-W(v)^(f£ ) -f¥). (4) 

The collision integral Eq. is only nonvanishing if v + and V- are nonzero simultaneously. It is localized at distances 
s d: The averages (is) decay at s > a away from the contact, being proportional to the solid angle at which the 
orifice is visible from the position point. Since is is constant within the orifice, t~ x J ds (is) ^ = 7^(0) where 7 ~ d/l 
is a contact-dependent geometric factor; we assume a constant 7 for simplicity. Finally, putting (f> = \i ~~ Xi we find 

= -7^(0K(0)(/£ 3 -4 13 ). (5) 
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FIG. 1: Spectral flow through the midgap states. For <f> > 0, a particle with p x > is captured at (j> = with an energy 
e+ = — |A| and then released at cf> = 2n with e+ = |A|; a particle having p x < is captured with e_ = |A| and released with 
e_ = — |A|. The dark dot (expanded in the insert) shows the impurity band ~ |A|y^y near the crossing point. 



SPECTRAL FLOW IN A BALLISTIC CONTACT 

In the collisionless limit 0], 



M0) = 7rV|A| 2 -e 2 <5(e-e±). (6) 

Since de±/d(f> — ±||A| sin(0/2) we can write Eq. (J5J as u±(0)df^'/dt — where the total time derivative of the 
distribution function is taken along the spectrum Eq. Q). The total time derivative vanishes: the distribution is 
constant in time and thus it is independent of energy. To find it we note that, as <j) varies from to 2tt, the distribution 
of particles with p x > remains the same as it was at t = 0, <p = for an energy e = — |A|, i.e., /W = — f A where 
/a = tanh (|A|/2T). For particles with p x < 0, the distribution remains as it was at t = 0, = 0, and e = |A|, i.e., 
/W = /a- If 4> decreases from 2-k to 0, the signs are opposite, 

= T/Asign(0). (7) 

Excitations are accelerated by the applied voltage undergoing multiple Andreev reflections 0,0,0 before their energies 
grow by 2|A| to let them escape from the potential well. 

The same process can also be viewed differently. Under a small voltage bias eV -C |A|, quasiparticles are accelerated 
along the spectrum Eq. (JJJ with = 2eVt. As the phase varies from to 27r, the quasiparticle energy changes from 
=F| A| to ±|A|. A quasiparticle that has been captured below the gap edge emerges above the gap after one period of 
Josephson oscillations (see Fig. QJ. Another particle gets captured and then released during the next period, and so 
on ■ This is known as the spectral flow. It plays a crucial part in vortex dynamics [l3l llif . In a context of weak 
links, it was first discussed in Ref. [lTj. In the absence of collisions, the distribution of excitations localized near the 
contact Eq. Q is determined by distribution of excitations at the gap edges, which are delocalized and are thus nearly 
in equilibrium due to the relaxation in the bulk. However, emerging quasiparticles have a "wrong" distribution with 
respect to the "native" ones; they relax producing dissipation. 

The current through the contact with the Sharvin conductance R^ 1 — e 2 NpVFS/2 is 



1 f de r 



r(l) ,(1) 

v+J+ -v-jl' 



(8) 



eR N J 2 

where S ~ a 2 is the area of the contact. With Eq. Q the current due to Andreev states is 

|e| J R J v/ = 7r|A||sin(0/2)|/ A sign(y). (9) 

Excitations with energies |e| > |A| give a much smaller current, which can be neglected. Equation agrees with the 
results ofRefs. 0,0,0. 



SPECTRAL FLOW AND IMPURITIES 

A small amount of impurities in the area of the contact lifts the degeneracy of Eq. JIJ and modifies the spectrum in 
a way that an impurity band appears near the point <f> = it where the two spectral branches cross 0] . In the presence 
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of impurities the DOS is i/ £ ±(0) = {E T /2^) ^(Aj/E+E-) - 1 where B± = e/\A\ ± cos(0/2). Both and are 
nonzero near the lines i?± = within the impurity band < E + E_ < 47 shown in Fig. ^ To simplify the partial 
differential equation © we note that, in the region far from the crossing point, |e| >• |A|^/7, \<f> — tt| » ^7, the ratio 
V-/v + = E + /E_ is of the order of -7I A| 2 /e 2 <C 1 near the line E + = 0, and it is v~/v + ~ e 2 /7|A| 2 ^ 1 near i?_ = 0. 
Therefore, far from the crossing point, \cf> — ir\ > (3 where y/7 <C /3 <C 7r, the DOS has the form of Eq. and the 

l.h.s. of Eq. 10 reduces to the full derivative df±^ /dt along the respective lines E± — 0. Since the collision integral 
also vanishes for \<j) — 7r| > (3, the distribution is constant in time far from the crossing point, in the same way as it 
was without impurities. 

Equation (JSJ) can be solved by perturbations in the limit of small voltages. The zero-order term corresponds to 
vanishing of the l.h.s. and results in the condition /i 1 ' = f^' that holds through the crossing point for a given energy. 
We find that the distribution functions at the four lines in Fig. entering the crossing region satisfy f± (n + (3) — 

/4 (tt ~ P) — i/A- To the zero approximation, particles emerging in the continuum with an energy e = ±|A| have 
the distribution ±/a that corresponds to the local equilibrium at the gap edges 

/i cq) = Tsign[cos(0/2)] / A . (10) 

In this approximation, no relaxation in the continuum occurs. The deviation f\± from the local equilibrium is 
proportional to the time derivative, i.e., to the voltage. We write f±~* = /± cq ' ) + /i±. The l.h.s. of Eq. JSJ) can 
be estimated as v eV (/ A /e*) where e* ~ |A|«/7 is the characteristic scale of energy near the crossing point. Using 
v ~ 1/^/7 we find fi± ~ =p (eV/|A|7) /a- The deviation f\± is thus linear in V for low voltages. 

In general, f±^ can be found by matching the four constants /I 1 (7r±/3) while integrating Eq. JSJ) across <j> = ir. We 
note that for a contact with a barrier of reflectivity R = 1 — D, the excitation spectrum is given by Eq. (J2J) with a gap 
\/i?|A| that separates states with positive energies from those with negative energies. As a result, the distribution 
for excitations with e > is decoupled from that for e < unless the Zener tunneling mixes them 0, . On the 
contrary, in macroscopic contacts with impurities, the states near <fr = ir form a continuum where the regions e > 
and e < are connected along the lines of a nonzero DOS passing through <fi = n, see Fig. Q] According to Eq. JSJ, 
the distribution of excitations is thus continuous within the entire region of nonzero DOS. Therefore, the distributions 
on all four lines entering the crossing region in Fig. ^ are coupled to each other. 

Since we are interested the asymptotic values of f± for <j> = ir — (3 and <f> = ir + f3, we neglect the fine structure of v 
and approximate it by Eq. © . The collision integral is then nonzero only at the crossing point of e+ and e_ defined 
by Eq. (JIJ . Equation JSJ) yields 



= T^\A\ 2 TS(e-e T ). (11) 



dt 

r is a combination /l 1 ' (e, <ft) — /_ (e, <j)) taken in the vicinity of e = 0, <fi = ir, which couples the four values f±^ (tt±{3). 
We shall see that r varies slowly near <j> = tt thus we approximate 
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(tt + P) + /I 1 ) {it -13)- (tt + 13) - f W (tt - 0) 



We integrate Eq. i|ll|) over time along the dependence E + = for particles with p x > or along £L = for p x < 
and obtain 

/i 1) ( 7 r + /3)-/£ 1) (7r-/3) = Tar (12) 
where a = Tr\A\-f/2eV. Adding the two Eqs. (|12(l we find that the quantity 

A(/3) = 4 1) (7r + /3)-/i 1 V-/?) (13) 

satisfies A(j3) = A(—(3) = A, i.e., it is a slowly varying function near (3 = 0. This proves that T — [A((3) + A(—f3)]/2 
can indeed be treated as a constant, T = A. 

Consider first the case <fi = 2eV > 0. The initial condition at t — 0, — is the equilibrium distribution at 
an energy e = =p| A|, i.e., — =f/a- Since the distribution is constant outside the crossing region, we have also 
f { ±\n - (3) = =F/ A . To find the functions /P(tt + (3) we use Eqs. (JHEJ whence A = -2/ A / (1 + a) and 

/i 1) (7r + /3) = T [(!-«)/(! + a)] /a- 
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For weak scattering a -C 1, the distribution remains unchanged as compared to its initial value, f± '(-tt + (3) = 

f±\ft — (3) — =f/a- Quasiparticles relax towards the local equilibrium ±/a as they emerge in the continuum from 
inside the gap at an energy ±|A|. If scattering is strong a 3> 1, the distribution changes considerably at the crossing 
point <j) — tt due to exchange between particles moving from the opposite sides of the contact. Emerging particles 
with an energy e = ±|A| have a distribution close to the "local equilibrium" Eq. (|10|) such that no relaxation occurs. 
This is not the true equilibrium for |e| < |A|: the latter can only be established when the spectral flow rate eV is 
much smaller than the inelastic rate. The local equilibrium is achieved through a very effective branch mixing due to 
the momentum exchange on impurities when particles pass through the crossing region of the spectral branches. The 
deviation from the local equilibrium is zero as long as <fi < tt while it appears for <f> > tt: 

fx± = /£ 13 (tt + /?) - [±/a] = T [2/ (1 + a)} /a ■ (14) 

This agrees with our earlier perturbative estimate for the low- voltage limit. 

Calculating the current for <f> < tt from Eq. (jSJ) we find I — I s where the "supercurrent" is defined through the 
"equilibrium distribution" Eq. i(TT)|l : 

eR N I s = tt\A\ sin(0/2)/ A sign [cos(0/2)] . 

This expression differs from the well-known result for a steady-state supercurrent 3]. The reason is that Eq. (|l(Jfl is 
not the true equilibrium, see the discussion above. For <fi > tt we have / = /,+/„; the "normal" current contains the 
deviation from equilibrium Eq. (|14[). Considering also the case <p = 2eT/ < we find for the normal current 

\e\R N I n = 2tt\A\ (1 + Hr 1 |sin(0/2)| / A sign(V). 

As in Ref. 0, it exists only when tt < <j> < 2tt for 4> > or when < cj) < tt for </> < 0. 
For a voltage-biased contact, the supercurrent vanishes after averaging over time: 

\e\R N l = 2|A| (1 + H)- 1 / A sign(V). (15) 

For high voltages, \a\ ~ |A|7/|eV| < 1, the MAR result Eq. © is reproduced. For low voltages |A|7 > \eV\ > t" 1 , 
the current is linear, / = GV where G = (4/7T7i?Ar)/A- This conductance is much larger than the contribution ~ i?^ 1 
from the states with |e| > |A|. To our knowledge, this low-voltage region has not been studied experimentally though 
the contacts with suitable parameters are now under intensive investigation fli| . 



DISCUSSION 

The linearity of the I-V curve for low voltages disagrees with a square-root behavior predicted in j^| using the 
model based on Zener transitions within independent channels averaged over the universal distribution of reflection 
coefficients valid for diffusive conductors. This distribution cannot be applied directly for nearly ballistic contacts. 
The proper distribution for contacts with a small disorder was found to be P{R) = (27T7)- 1 % /7 Tr l?/v A R. Its 
denominator does also contain y/R like the universal distribution, which, at the first glance, should have also resulted 
in a square-root current-voltage dependence. However, the kinetic-equation approach applicable for a multichannel 
contact leads to a different result. The reason is that, in a macroscopic contact, the quantum channels arc very 
close to each other in energies with the level spacings inversely proportional to the number of channels N c h_. Thus, 
in addition to the transitions within each channel (which are the only ones considered in the independent-channel 
model), time dependent processes excite even more effectively transitions between different channels. As a result, the 
distribution function is spread over the entire energy interval of nonzero DOS or over the range R ~ 7 of reflection 
coefficients. The role of each particular channel within the range SR should thus be weighted with SR/^f. If we now 
average the low-voltage Zener processes having the rate oc exp(—\A\R/eV) with P(R) we would obtain the essential 
reflection coefficients R ~ eVY|A|. Thus the fraction of essential channels contributing to the current in multichannel 
contacts is decreased by a factor eV/|A|7 as compared to the independent channel model. This yields a much smaller 
current compared to what is predicted by the kinetic-equation approach. Therefore, the Zener processes can be 
neglected in favor of the momentum relaxation on impurities. To summarize, the model based on Zener transitions 
within independent channels docs not seem to be suitable for kinetic processes in multichannel contacts. We believe, 
however, that calculations using the Landauer approach with the full account of interchannel transitions followed by 
the appropriate averaging would be equivalent to the kinetic equation technique. 
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